A variational approach to necklaces formation in polyelectrolytes 
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By means of a variational approach we study the conditions 
under which a polyelectrolyte in a bad solvent will undergo a 
transition from a rod-like structure to a "necklace" structure 
in which the chain collapses into a series of globules joined by 
stretched chain segments. 



I. INTRODUCTION 

It is well established that in a salt-free environment 
and under good or 0-solvent conditions, the electrostatic 
interaction between the monomers of a single polyelec- 
trolyte chain induces a rod-like conformation. In the 
pressence of of a small amount of salt, if the ions are ho- 
mogeneously distributed such that their interaction with 
the chain is limited to screen the elctrostatic interactions, 
a rod-like conformation is also expected. Even in this 
simple model, the introduction of a a secondary inter- 
action can leai|to interesting structural modifications. 
Dobrynin et aZ.EJ have presented a scaling argument and 
computing evidence that by turning on an attractive in- 
teraction between the monomers, the standard rod-like 
conformatiorJj, is abandon in favor of a "necklace" one. 
The necklace structure consists of highly stretched seg- 
ments alternating with collapsed globules, the "beads" 
of the necklace. Another mean field calculation by Bal- 
azs et ala also produce similar results. In this paper 
we further explore this possibility, by means of a varia- 
tional approach. The techniques we use here are adapted 
to the case of long (effectively infinite) weakly charged 
polymers. 

The necklace structure was first suggested by Kantor 
and KardarQ in their study of polyampholytes. This was 
argued in|-|the basis of a similarity with the Rayleigh 
instabilityO that appears in charged fluids. Namely, 
a large charged droplet of fluid will split into smaller 
droplets to put the charges far away from each other and 
the splitting is only stopped by the increasing cost of sur- 
face energy. The droplets are kept together by cohesional 
molecular forces, which in the polymer case are replaced 
by the bonds between monomers. 



II. VARIATIONAL APPROACH TO 
POLYELECTROLYTE CONFORMATION 



this and other techniques has been reviewed of Barrat 
and JoanyQ). p 

The basic variational principleQ states that given the 
Hamiltonian of the system if, and a family of simpler 
trial Hamiltonians Ht characterized by a set of param- 
eters a,b, . . ., the free energy of the system satisfies the 
inequality 



F<Ft + {H- Ht)t 



(1) 



where Ft is the free energy of the trial Hamiltonian, and 
the braces indicate averages with respect to the statisti- 
cal ensemble associated with the trial Hamiltonian. From 
among the family of trial Hamiltonians one can then 
choose the one that gives the lowest upper bound for the 
true free energy. Approximations to the expectation val- 
ues of observables of the system can be then taken to be 
those calculated by using the optimal trial Hamiltonian. 

The usual problem with this method is that while the 
bound on the free energy is rigorously true, it might be 
useless and misleading if the family of Hamiltonians is 
not close enough (as operators in some space) to the true 
Hamiltonian. Solving this problem requires the use of 
general enough trial Hamiltonians with which practical 
calculations can still be done and that from a heuristic 
point of view one can expect them to capture the essential 
features of the system. Further, a successful application 
of this method can be very instructive if the parameters 
that describe the trial Hamiltonians have a direct physi- 
cal interpretation. 

The model Hamiltonian of the polyelectrolyte chain 
that we consider is the well established (continuum ver- 
sion) form, 
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dn'V{r{n) - r(n')) 
(2) 



with the usual entropic term and with the interaction 
potential between segments V modeled as having two 
parts: a hard-core interaction 



Vhc = i^o<5(r) 



(3) 



that we will consider attractive (i.e. vq < 0), wich corre- 
sponds to the immersion of the ploymer in a bad solvent, 
and a screened electrostatic part 



The variational approach has been often used in the 
investigation of the structure of polyelectrolyte chains. 
This technique was first introduced in this context by de 
Gennes et ala. (Relevant work in this area by nreans of 



Ve = fel- 



er 



(4) 



where / is the valence per monomer, eo the electronic 
charge, and e the permittivity of the medium. We will 



consider only the case in which the charges are quenched 
in the polymer backbone. 

To observe the necklace conformation, the trial Hamil- 
tonians must explicitly include the possibility of uniform 
and non-uniform stretching. We choose the following ver- 
sion of the stretched chain: 
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(5) 



This trial Hamiltonian incorporates Gaussian fluctua- 
tions around a background stretched conformation Tq 
with a renormalized Khun length b'. In the basic 
stretched chain model, the conformation has uniform 
stretching but we will consider strechings with further 
spatial dependence. 

More precisely, we will consider background conforma- 
tions of the form 



ro(n) = ppnu + e cos{qn)u 



(6) 



corresponding to a uniform unidirectional stretching 
modified by small modulations of the stretching along the 
same direction. This modulation in the average stretch- 
ing rp creates zones of accumulation and depletion of 
monomers along the basic rod-like conformation. We will 
show that for a range of values of the basic parameters 
of the problem, the minimum free energy with a uni- 
form stretching is larger than if we include modulations 
at certain wavelengths; i.e., we will not find the optimal 
stretching but prove that the rod conformation is unsta- 
ble against infinitesimal modulation. The wavelength of 
the most unstable mode thus predicts the spacing of the 
beads. (See figure 1). 



Figure 1 . 
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FIG. 1. Typical configurations of the polyelectrolyte 
chain: (a) uniform stretching and (b) under modulated 
stretching with wavelenght 2-nb/q. 

Bratko and DawsonO have considered the same model 
Hamiltonian and a similar variational principle that 
uses a general Gaussian matrix of correlation between 
monomers and a constrain on the allowed conformations 
that keeps the bond-length fixed. This last constraint 
has also been incorporated in a variational approach by 



Ha and Thirumalaia in their disccusion of the rigidity 
of a polyelectrolye and, as discussed below, regularizes 
the behavior of the model in the strongly stretched and 
collapsed limits. In the previous studyEl of the polyelec- 
trolyte in bad solvent a persistent stretching was not con- 
sidered {i.e. dvp/dn is set to 0) and the stretching of the 
chian was only observed through the large amplitudes of 
the fluctuations (r(ri)r(n')). Instead, our approach takes 
the average stretched conformation rg as the basic quan- 
tity describing the system. 

The possibilities of uniform and non-uniform stretclt 
ing have been properly considered by Jonsson et oLEj 
when looking for the corrections to the stretching arising 
from the finite size of the polyelectrolyte. Their work, 
however, has not included consideration of further inter- 
actions beyond the electrostatic one. 



III. BASIC ROD-LIKE CONFORMATIONS 

We start by looking at the minimum of the varia- 
tional energy that can be obtained by considering only 
trial Hamiltonians with uniform background stretchings. 
Thus, we write 



P{n) = Po 



(7) 



and the expressions in (|I|) are calculated as follows. First, 
the free energy for each of the springs associated to 
each pair of neighboring monomers is, up to a constant, 
— SfesTlnfe'. Next, the average of the stretching en- 
ergy per monomer(i.e. per unit length of the continu- 
ous chain) of the original chain in this ensemble has the 
forced background contribution {3/2)kBP^ and the basic 
random walk part {i/2){h' /hfkBT. 

The probability distribution for the position of two 
chain segments separated by n monomer units G(r(s -|- 
n),r(s)) — G{r{n) — r(0)), is just a Gaussian function 
with mean 



(r(n) - r(0)) = ppnu 



and variance 



((r(n) - fo - pnuf) = b" 



(8) 



(9) 



Therefore, for an interaction potential V{r), the averaged 
interaction energy per unit length is 

{V) = 2 f ds f drG{r{s) - r(0))V^(|r(s) - r(0)|) (10) 

or using the potential's Fourier transform and the explicit 
form of the correlation between monomers: 

/•OO 

{V)^2 ds 
Jq 

■exp[ \-ik-upos]V{k). (11) 



(27r)3/2 
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Integrals of this type that diverge on the s -^ hmit, 
are regularized by the short distance cut-ofF Smm — 1 
(namely, one bond distance). 

The final result of the evaluation of the averages is (in 
units of ksT), 



F = -i\nh'+'-{h'^ + pl 



I) 



In 1 + -n^ 

Po V b'^K 




(12) 



In this expression and in the rest of the paper, all dis- 
tances are measured in Kuhn length units. We have also 
introduced the Bjerrum number zb ~ f^e^ /kBTeb, and 
a dimensionless interaction pre-factor i^ = lygb^ /ksT. 

Minimization of expression (^2|) with respect to po, and 
6', gives the basic stretched conformation as a function of 
the basic parameters, k, zb and v. In general, there exists 
two minimums for this expression: one in the stretched 
region pQ ^ and another in which the chain is collapsed 
(or Gaussian) po ~ 0- The present method and approxi- 
mations, however, is suited to study only stretched con- 
formations and cannot be reliably used to decide between 
the two states. For example, in our proposed variational 
family of Hamiltonians the renormalized Kuhn length b' 
is intended to take values in the neighborhood of 1 (i.e. 
b' « b). Thus, to study the globule-like collapsed chain 
it is necessary to modify the present variational princi- 
ple and introduce a slightly different model Hamiltonian 
(by, for example, introducing a bond length constraint) 
that prevents the unphysical collapse &' ^ 0. To a large 
extent results concerning this collapsed|-|regime can be 
found in the work of Bratko and Dowsodl. The possibil- 
ity that a collapsed state might have a lower energy than 
the stretched state leads to modifications of our results 
that will be pointed out in the rest of the paper. 

In finding the minimum of the evaluation of the av- 
erages, eq. (|l^), we observe that for large values of the 
attractive interaction there is no local minimum with fi- 
nite stretching. In figure 2 we show a numerical example 
of the behavior of the stretching and the renormalized 
bond length, as functions of i^ and zb- According to the 
previous discussion, the real stretched region in the v-zb 
diagram is smaller than the one shown since the collapsed 
minimum will have a lower energy when the stretched 
minimum becomes unstable. In the examples provided 
we use small values of the Bjerrum number, to keep the 
stretching consistent with our assumptions of Gaussian 
springs. Stretchings of order one should be treated with 
a suitable representation of the constraint of fixed bond 
length. Further, in the units we use, we look at regions 
of parameter space that corresopond to fairly large hard 
core interactions but that still represent experimentally 
realizable conditions. 



Figure 2a. 
Uniform Stretcliing. 




Figure 2b. 
Renormalized Khun lenght. 




FIG. 2. The ~v~zb diagram for k = 10"'', where the thick 
lines mark the boundary between the collapsed and rod-like 
states. Within the rod-like region it is shown (a) lines with 
constant values of the uniform stretching po and (b) lines with 
constant values of the renormalized Kuhn lenght. 

Finally, we point out that the criteria for the validly of 
these results in the case of finite size chains is as follows. 
The induced end to end stretching p^N must be larger 
than the natural chain size N^'"^ . Thus, our conclusion 
will be relevant for systems satisfying the condition 



Po > N''/^ 



(13) 



IV. 



STABILITY OF ROD-LIKE 
CONFORMATIONS 



We will now determine the stability of the solutions 
found for uniform stretching. For this, we introduce a 
small perturbation on the background stretching, with a 
well defined wavenumber q measured along the backbone: 



p(s) = pos + PqCOs{qs). 



(14) 



ith pq being a small amplitude. The variational free en- 
ergy per monomer associated with this stretching can be 
expanded in a power series in the amplitude pq 



f = .fo 



Plfq 



(15) 



Obviously, the uniformly stretched conformation does 
not have the minimum energy if the fq is negative. The 



true minimum can be found by minimization over all 
possible stretchings, but we will not attempt to solve 
this problem here. The results we obtain, however, show 
that when the instability occurs, it does so in a partic- 
ular range of wavelengths, indicating that the minimum 
energy is obtained in a regular necklace conformation, 
rather than in a randomly spaced one. 

The linear stability analysis against small perturba- 
tions can also be carried out for perturbations that 
are transversal to the chain. These perturbations do 
not lead to instabilities but their study is important in 
understaading the persistence length of palyelectrolyte 
chainalirO'EJ and will be studied elsewheral3. 

The contribution to the free energy correction from 
the entropic Gaussian part of the model Hamiltonian is 
simply (3/4)(7^/f)^. For the potential interaction part we 
write the contribution of a pair of monomers with labels 
t' and t {t' > t) as 



/ 



dk 



(277)3/2 



exp[ 



Q{t' - t)fc2 
5/2 



exp[zk • (p(t') - pW)u]T^(k). 
(16) 



that occurs for sufficiently large values of the attractive 
interaction disscused before, the rod-like conformation 
region in which the electrostatic interaction is dominant, 
and the more interesting part for our purposes, a narrow 
strip between them which is still rodlike but has under- 
gone a finite wave-length modulation. In this last region, 
we have identified the points for which the value of the 
critical wavenumber q is the same. We also note that the 
instability region for necklaces starts always at a finite 
value of the Bjerrum number, since it is necessary for the 
system to have enough electrostatic energy to elongate 
the chain. 



Figure 3. 
Instabifity Region. 



We expand the exponential factor into 




q* 

a 0.025 
b 0.045 
c 0.050 



0.004 0.006 



exp [ik • {p{t') — p(t))u] = exp[ik • po{t' — t)u\ 
(1 -I- ik • upq{cos{qt') — cos{qt)) 
1 
2^1^ 



-^plk ■ u^{cos{qt') - coa{qt)f + 



(17) 



This expansion is substituted into (16), a new variable 
s = t' — t is introduced and the integration over t elim- 
inates all oscillatory terms not depending on s only, to 
obtain the final result, 



2 



V,P\ 



1 



ds 



iP\l, -J (2^)3 



dk 



(k • u)2 [-1 + cos(gs)] exp[ik • pou]V"(k). (18) 



The criterion for instability is then: 

U = J?" + n < 0. 



(19) 



Concrete evaluation of eq. (jlj) can be carried out for 
specific values of z^, k and u. The criteria is met for 
a bounded region of q values. Among the q values for 
which the system is unstable there will be one value q for 
which the determinant function fq is most negative. We 
identify this wave-number with the inverse of the most 
likely spacing between the beads of the necklace. 

It is easier to discuss the results keeping k constant. 
For a given value of k we can compute the region of sta- 
bility in the zb-v plane and we present one such example 
in figure 3. On changign k, the general morphology of 
the diagram is still the same as in figure 3 and our re- 
sults are consistent with those obtained by Dobrynin et 
alh. There are three main regions. The collapsed state 



FIG. 3. Region of neckalce instability for k — 10~* in the 
diagram ~v — zb- Point a marks the limit of the necklace 
region, and lines b and c mark points of equal values of the 
most unstable wave-numbers q* . 

Different values of k produce the same general re- 
sults in that we always reach, first, a point of instabil- 
ity against collapse, and for sufficiently large charge, and 
higher values of the attractive interaction, a necklace pat- 
tern. The main effect of a decreasing k is to make the 
electrostatic interaction stronger, thus pushing the re- 
quired values of v for each instability ever higher. 

We add further detail to our example to illustrates 
another feature of the system. We have computed for 
K = 10~^ and ZB = .001 the region of instability and the 
region of unstable modes. The results appear in figure 4. 
There, we confirm the different nature of the collapsed 
and the necklace state. The region on q space which is 
unstable is always finite. To obtain a smooth transition 
to the collapsed state we would need this region to grow 
and touch large wavelength modes {q « 0) but this is not 
the case. 



Figure 4. 
q-values for instability 
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FIG. 4. Necklace unstable modes q for k = 10^ and 
zb = 0.01. The dashed line marks the most unstable mode 
for each value of u. 



V. CONCLUSIONS 

We confirm necklace formation in polyelectrolytes us- 
ing the variational principle technique. In this calcula- 
tional scheme we are limited to consider the regime in 
which the entropic free energy of the chain is still com- 
parable with the energy scale imposed by the charge in- 
teraction, which leads to weak stretching. This calcula- 
tion gives us not only insight on the onset region of the 
phenomena, but it also provides a concrete parameter 
space window for the experimental search of the tran- 
sition. The present calculation can be refined in many 
ways. As we have mentioned before, it is important to 
consider the effects of a fixed bond length, and secondly, 
we can complicate the Gaussian correlations in the trial 
Hamiltonian so as to consider non-neighboring monomers 
correlations. 

In the future it will be very important to determine 
the exact location of the transition between the collapsed 
and extended conformations. Since the region where the 
necklace formation is observed lies between these two 
states, it might happen that the collapsed state acquires 
the lowest energy before the onset of the finite wave- 
lenght instability. In that case the necklace cannot be 
observed. On the other hand, this scenario should not 
be true in general since our method can also underesti- 
mate the point at which the necklace transition occurs. 
Again, the physical onset of the necklace preceeds the fi- 
nite wavelenght instability of the rod-like conformation. 

Finally , from the methodological point of view, we 
believe that the uses of the variational approach have not 
been exhausted. It provides means to test the different 
structures arising in single and multiple chains systems. 
For example, ona can study the case in which the charges 
are not quenched aloiM the chain, a problem considered 
by Higgs and Joannyll^ and recently jextended to study 
interchain interactions by Liu and Hall3. 
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